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ABSTRACT 


Images of an airborne, scanning, radiometer operating at a 
frequency of 98 GHz, have been analyzed. The nun-wave images were 
obtained in 1985/1986 using the JPL nun-wave imaging sensor. The 
goal of this study was to enhance the information content of 
these images and make their interpretation easier for human 
analysis. In this report, a visual interpretative approach was 
used for information extraction from the images. This included 
application of nonlinear transform techniques for noise reduction 
and for color, contrast and edge enhancement. Results of the 
techniques on selected mm-wave images are shown. 
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CHAPTER I. INTRODUCTION 


Passive millimeter-wave sensors are very useful for 
battlefield surveillance because they provide information during 
foggy and cloudy weather when visual and IR sensors are unusable. 
This sensor can be used from an aircraft, helicopter, or from a 
ground based platform. One of the main advantages of the passive 
mm-wave imager, is that the interpretation of these images by 
human operators is similar to the interpretation of visual and IR 
images because of the high correlation between the intensities of 
adjacent points. This is in contrast to radar images which are 
complicated by speckle, polarization and edge backscatter 
effects. 

Passive sensors have the advantage that the observations can 
be made without transmitting radiation. This reduces the 
possibility of being detected, destroyed, jammed or fooled as can 
happen with active radar sensors. However, radar sensors can 
provide high- resolution range and Doppler information which aids 
in target identification. Thus, a good surveillance system would 
use a passive sensor to survey the overall scene and then use the 
radar sensor in short bursts to help identify targets. A 
comparison between these passive and active sensors is summarized 
in Table 1-1. 

One of the limitations of the passive mm-wave sensor is that 
it does not have the fine spatial resolution of visual and IR 
imageries. To help in the interpretation of this lower spatial 
resolution data, this study on mm-wave imaging data acquired by 
JPL in 1985 and 1986, was done to develop and apply special image 
processing techniques for the mm-wave images. There are two 
approaches which can be used for image information extraction. 

One approach, which makes use of the computer, is known as 
quantitative analysis or classification. The other approach, 
based on a visual inspection of an image by a human analyst, is 
known as photo or image interpretation. The properties of these 
approaches are summarized in Table 1-2 t 1 ]. Based on our 
experience with the mm-wave images, it was concluded that the 
photo interpretation approach is best for the evaluation of 
passive mm-wave images for battlefield surveillance. The shapes 
of important large scale features like roads, coastlines, rivers 
and other lineaments, which can be derived from mm-wave images as 
high contrast edges, are easier to interpret by a human analyst 
than by a computer. On the other hand the quantitative analysis 
approach can be considered complementary to aid in future auto- 
matic data analysis. 
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Table 1-1. Surveillance Sensor Comparison 


Parameter Passive - Radiometer Active - Radar 


Spatial Large: 5-10 m 

Resolution 

Identification Detection - Limited 

classification 

Range/Velocity Small: 1 km wide swath 

Range: <2 km 


Image Easy to understand. 

Interpretation similar to visual - 

changes in emissivity. 
Correlates with visual 
maps . 

Not sensitive to aspect 
angle and polarization. 

Real-time analysis 
possible. 

Computer Simple image processing 

Processing techniques can be used 

to assist observer. 


Data Rate Low: -100 kHz 

Security Undetectable - no 

emission. 

Difficult to jam - need 
kWatts. Difficult to 
spook, hard to 
camouflage targets. 


Equipment Relatively simple and 

inexpensive. 


Small: 1 m 


Detection and 
Recognition 

Large: -5 km wide 

swath 

Range: -10 km 

Harder to understand 
because of speckle, 
aspect angle, depres- 
sion angle and 
polarization. 

Major returns from 
edges and corner 
reflectors . 


Large data set to 
analyze. 

Requires sophisticated 
computer processing. 

High: -1 MHz 

Detectable RF emis- 
sion, can be locked 
on and destroyed. 

Can be jammed and 
spooked . 

Can be fooled by 
camouflage, corner 
reflectors and 
transponders . 

More complicated and 
costly. 
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Table 1-2. A Comparison of Approaches of the 
Interpretation of Digital Image Data 


Photointerpretation Quantitative Analysis 

(by human analyst) (by computer) 


Occurs on a large scale relative 
to each pixel 

Inaccurate area estimates 


Limited multispectral analysis 


Assimilate only a limited number 
of distinct brightness loads 
(«16 levels) 


Shape determination is easy 


Spatial information is easy to 
use in a qualitative sense 


Occurs at individual pixel 
level 

Accurate area estimates 
possible 

Multispectral (multi- 
dimensional) analysis 

Uses quantitatively all 
available brightness 
levels in all features 
(e.g. 256) 

Shape determination 
involves complex software 
decisions 

Limited techniques 
available for making use 
of spatial data 


Image processing techniques which were found to be most 
useful, especially in real-time analysis, were the choice of the 
false color map, modification of the brightness temperature 
histogram for contrast enhancement, smoothing for noise 
reduction, and edge enhancement. These techniques are discussed 
in the following sections and examples using selected nun-wave 
images are shown. 
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CHAPTER II. BACKGROUND 

When a mm-wave sensor observes a ground scene from 
above, as shown in Fig. 2-1, the received signal power or 
brightness temperature consists of thermal emission from 
objects in the antenna beam, reflected sky emission from 
these objects, and atmospheric emission from the atmosphere 
below the sensor: 

Tree = T 0 bj + T r ef sky + Tb c ]^gnd + T a tm 



Figure 2-1. Millimeter-wave radiometer observing a 
ground scene from above 


The received brightness temperature from objects in the 
antenna beam is equal to the product of its physical 
temperature and its surface emissivity, multiplied by the 
ratio of the object's area to the antenna's beam area and 
reduced by the atmospheric attenuation: 

T ob j = e*T*n/L a (K), (2) 
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where : 


e emissivity of the object and is a function of 

the material and its surface roughness 
(0 < e < 1) , 

T physical temperature of the object (K) , 

L a atmospheric absorption factor below the 

radiometer, i.e., = 10 E (L a t m *R/10000 ) > 

L atm atmospheric attenuation (dB/km) below the 

radiometer, 

R range to object (m) , 

n ratio of object's area to antenna beam area 

(0 < n < 1), i.e., A/ ( P 1/ 4) ( R-t an® b ) • 

A object's area (m2) 

e b antenna beamwidth (deg) 

The reflected sky emission temperature is equal to the 
product of the radiometric sky temperature and the quantity 
(1 - € ) and is also multiplied by the ratio of the object's 
area to the antenna ' s beam area and reduced by the atmo- 
spheric absorption: 

T r ef sky » (1-e) *T sky *n/L a (K) , (3) 

where : 

T skv radiometic sky temperature 

= T atm *(l - 1/L Z ) + T cb /L z (K) 

T atm mean atmospheric temperature ( =: 280 K) , 

L z total zenith absorption factor, 

T cb cosmic background radiation temperature 

( - 2.7 K) . 

The received signal temperature from the background 
surrounding the object is given by: 

Tbckgnd = ( e g* T g + i 1 ~ c g) * T sky)/ L a > * ~ *) ( K ) ( 4 ) 

where is the emissivity of the background and Tg is the 
physical temperature of the background. The atmospheric 
emission below the radiometer is: 


5 



Tatmos = T a *d - 1/L a ) (*> » (5) 

where is the physical temperature of the atmosphere below 
the radiometer. 

MM-wave scenes can show large contrasts with reflecting 
objects, such as metal vehicles, buildings, or wet areas ap- 
pearing much colder than the surrounding rough absorbing 
background. Thus, ref lecting objects will "stand out" from 
the surrounding background with a large temperature 
difference or signal-to-noise ratio. 

The system configuration used with this mm-wave imager 
is shown in Fig. 2-2, and a simplified block diagram is shown 
in Fig. 2-3. A complete technical description of the mm-wave 
imaging system was published in Ref. (2) and the characteris- 
tics are summarized in Table 2-1. In the mm-wave imaging 
sensor, a mechanical scanning system scans the 0.5 s antenna 
beam ±22° across the flight track. This beam is stabilized 
to correct for aircraft movement due to turbulence and uses a 
butterfly scanning pattern to produce a linear raster scan 
pattern on the ground as shown in Fig. 2-2. This feature 
allows the mm-wave image to be displayed in real time. The 
operation frequency of 98 GHz was chosen to give the smallest 
spatial resolution with the best poor weather performance. 

An altitude of 750 m was chosen as a typical value for 
battlefield surveillance and provides a ground resolution of 
7 m. If larger spatial resolutions can be used, the 
operating altitude could be increased. 
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Figure 2-3. Simplified block diagram of the mm-wave 
imaging sensor 


Table 2-1. JPL Millimeter-Wave Imaging Sensor 
Characteristics 


Aircraft Altitude 

750 m 

Aircraft Speed 

41 m/s (kts) 

Scan Width 

±22.5° or 620 m 

Angular Resolution 

0.5° 

Ground Resolution 

<7 x 7 m 

Radiometer Type 

Total Power, Solid-s 
Dual Polarized 

Center Frequency 

98 GHz 

Sensitivity 

<0.8 K 

Data Rate 

2800 samples per sec 

Data Display 

False color TV image 

Weight 

70 kg 

Power 

250 Watts 
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CHAPTER III. IMAGE PROCESSING TECHNIQUES 

There are numerous types of image processing techniques 
which have been developed and described in the literature (Refs. 
3, 4, 5, 6). The goal of this study was to select the best 
algorithms for the processing of the lower spatial resolution 
high contrast mm-wave imaging data. These algorithms had to aid 
in the data analysis and be fast enough to operate in a real-time 
environment. After considerable study, a program package was 
developed for processing the mm-wave data. A summary of the 
image processing programs developed during this study and a 
detailed list of the menus are in the Appendix. The description 
of the techniques used on selected images is presented below. 

The images in this report were selected because they provided 
good examples for the demonstration of these techniques. 

A. FALSE COLOR ENHANCEMENT 

The motivation for using color in image processing is the 
fact that the human eye can discern thousands of color shades and 
intensities. This is in sharp contrast with the eye's relatively 
poor performance with gray levels. Only about 16 shades of gray 
are detectable at any one point in an image by the average 
observer. Therefore, a color display has the capacity to convey 
much more information than an achromatic display and has the 
potential for improving the interpretability of images. 

A careful selection of the false color scale for the 
brightness temperature range is an important item in the 
interpretation of the mm image. Letting the observer adjust the 
color scale and temperature range in real time gives excellent 
results in this analysis. One way to think of this is that each 
observer has a matched filter in his brain to process 
information, and being able to change the color scale to match 
this filter, will maximize the information transfer 
between the observer and the computer. 

The color scale used in this report was generated 
empirically incorporating our former experience in data display. 
Normally colors can be represented^by^three dimensional vectors. 
If C represents a given color and R, G, B represent the unit 
vectors of the three fixed primaries, then the equation 

C = RR + GG + BB 

is used for the production of an additive mixture of the 
quantities R, G, B. In our case the tristimulus values R, G, B 
were described by harmonic waves with varying amplitudes, 
frequencies and phases. The calculated and optimized transfer 
functions for R, G, B are presented in Figure 3-1 and the 
resultant color scale is shown in Figure 3-3. More sophisticated 
procedures which can also be used are described in Ref. 7. 


9 



COLOR INTENSITY 



Figure 3-1. Transfer functions for the color map shown 
in Fig. 3.4 vs. pixel color intensity 
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The image of the 1-5 freeway near Carlsbad, Ca. on a 
partly cloudy day is shown in Fig. 3-2. Fig. 3-2a is the 
visual image and Fig. 3 -2b is the original nun-wave data with the 
false color scale adjusted to provide a realistic looking 
picture. This image is repeated in Fig. 3-3b, with a change in 
the color scale and gives a completely different impression of 
this scene. Choosing the color scale to appear similar to visual 
images seems to be one of the best choices to give the observer 
the best information transfer. Another example of the change in 
color scale is shown in Figs. 3-4b and 3-5a, in an image of the 
Agua Dulce Airport. The more realistic looking scene in Fig. 

3-5a is easier to understand and see the individual objects. 

The key to selecting the best color scale is the ability to 
do this quickly and easily. In our current data reduction 
system, the color scale offset (the position of the color map on 
the brightness temperature scale) and the gain (the spread of 
colors vs brightness temperature) are quickly adjusted 
with a track ball to give the best image. Thus in a few 
minutes the best colors can be selected by each observer, and 
only minor changes are required for different images. 


B. SMOOTHING 


Smoothing algorithms reduce the random noise and other 
spurious effects that may be present in an image as a result of 
sampling, quantization or disturbances during image acquisition. 
These algorithms provide a better signal-to-noise ratio 
especially for distributed targets. However, the spatial 
resolution will be degraded and the edges are blurred. Three 
different smoothing algorithms were used on the original mm-wave 
ship image shown in Fig. 3-6a: a moving average, a median filter 

in the spatial domain, and a low pass filter in the spatial 
frequency domain. 


the 


The moving average filter (or neighborhood averaging) 
form: 

n-1 m-1 

S'(x, y) = ( 1/ (n*m) ) 2 2 S(x + k, y + 1) , (6) 

k=0 1=0 


has 


where S(x, y) is the input image and S'(x, y) is the output 
image. The moving average filter with a 3x3 window, (n = m = 3) 
was used on the original image and the result is shown in Fig. 
3-6b. As seen, the image is quite blurred especially at the 
edges . 

A median filter, in which the intensity of each pixel is 
replaced by the median of the intensities in a pre-defined 
neighborhood of that pixel, has the form: 

S'(x, y) = Median (S(x, y) , . . . , S(x + k,y +1)), (7) 
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Fiqure 3-2. Video and mm-wave images of a section of the 1-5 freeway south of 
Carlsbad, CA, taken on a rainy day. The imager was flown over a 
layer of broken clouds as seen in (a). The area pictured in the 
mm-wave image (b) is 0.4 by 2.4 km. The image in (c) is contrast 
enhanced and the image in (d) has been processed with the Sobel 
edge detection algorithm. 
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This image is the same as in Fig. 3-4 with a different color map. The original mm-wave 
data is in (a), the contrast enhanced image in (b) and the edge detected image in (c). 







ORIGINAL PAdk? 

COLOR PHOTOGRAPH 
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Figure 3-6. MM-Wave images of a large ship in the Long Beach, CA. harbor just south of Terminal Island. 

The original data is in (a), and the smoothed data in (b) - (d) as explained in the text. 






with 0 < k < n-1 and 0 < 1 < m-1. The result of this algorithm 
(also with n = m = 3) is shown in Fig. 3-5c. This image shows 
the reduction in noise and that there is less blurring than in 
the moving average image. 

The third approach to smoothing was to low pass filter the 
image in the spatial frequency domain. The original image in 
frequency space, S(u, v) , was convolved with a low pass filter 
function H(u, v) to reduce the high frequency noise: 

S'(u, v) = H (u, v) * S (u, v) . (8) 

The low pass filter function had the form: 

H(u, v) = 1/(1 + 0.4* (Sqrt (u 2 + v 2 )/0 o ) 2n ), (9) 

where D 0 and n are variables used to minimize the blurring 
effects. The results of this operation, with D 0 = 4 and 
n-1, are shown in Fig. 3-6d. 

A comparison of Figs. 3-6c and 3-6d shows only small 
differences. Thus it was concluded that the median filter, which 
is easier to implement for real time processing, is the best 
smoothing filter to use in the mm-wave image processing because 
of the better edge preserving properties. 

C. BRIGHTNESS TEMPERATURE HISTOGRAM MODIFICATION 

Techniques of image enhancement have the goal to improve the 
visual appearance of an image and to convert the image to a form 
more amenable to human and machine analysis. Unfortunately, at 
the present time there exists no unifying theory of image 
enhancement since it is difficult to define a criterion 
characterizing a "good quality" image. One possible reason lies 
in the multiplicity of application fields. In our case, a 
technique which was found very useful was to equalize the 
brightness temperature histogram equalization or contrast 
enhancement (Ref. 4). This method is described in Fig. 3-7. The 
grey value histogram of the original image, which may have many 
of its pixels in a narrow temperature range, is mapped into a 
histogram with a nearly uniform distribution across the entire 
temperature range. After this modification, the cumulative 
distribution would be nearly a straight line. This means that an 
initial grey level interval with a high population will cause a 
larger stretch of the output grey levels than an interval with a 
lower population. Examples of images using this technique are 
shown in Figs. 3-2c through 3-5c, Fig. 3-8c and Fig. 3-9c. 

This processing enhances larger areas while suppressing 
smaller objects. This resultant image is especially useful 
for enhancing distributed targets and identifying large scale 
geographic features such as types of vegetation, roads, and 
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HISTOGRAM MODIFICATION 


t t t . t 



BRIGHTNESS TEMPERATURE 
(COLOR SCALE) 


Figure 3-7. Diagram explaining the brightness temperature 
histogram, H(T) , modification for the contrast 
enhancement algorithm 
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CUMULATIVE DISTRIBUTION 



hills. An example of this is Fig. 3-8 which is an image of 
military vehicles in an open field at Camp Pendleton. . The 
small variations in the land area next to the road, Fig. 
3-8b, are clearly emphasized in the contrast enhanced image 
in Fig. 3 -8c. However, this technique suppresses small 
targets with only a few pixels. Thus, the vehicles in this 
image are difficult to detect. 

Another good example of this technique is with the 
cloudy freeway scene in Fig. 3-lc. In this image the 
contrast enhancement makes it easy to see the freeway, the 
adjacent roads and the large fields near the freeway. 
Identifying these large areas is useful in determining the 
exact scene location. In the Agua Dulce Airport image, Fig. 
3-4c, the changes in the terrain slope are easily scene. 
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D. EDGE DETECTION 

In nan-wave images, boundaries of objects appear as 
brightness temperature or intensity discontinuities. For the 
human observer, an object can often be recognized from only a 
crude outline. To enhance the small objects and edges with high 
brightness temperature contrasts in the mm-wave image, special 
edge enhancement algorithms are used. Two algorithms from the 
main class of gradient and template matching operators tried on 
this mm-wave data were the Laplace and the Sobel operators. The 
formulas for these operators are given in Table 3-1. These 
operators reduce the temperature of all the pixels, except those 
on the border of a temperature discontinuity. The temperature 
increase at the object is proportional to the degree of the 
transition. Thus cold spots turn into hot doughnuts. The 
advantage of this algorithm is that it enhances small details 
that may not have been obvious on a uniform background. 


Table 3-1. Laplace and Sobel Operators 

LAPLACE OPERATOR 



A2 


A8 

AO 

A4 


A6 



G(X,Y) = 4 * AO — (A2 + A4 + A6 + A8) 


SOBEL OPERATOR 


A1 

A2 

A3 

A8 

G 

A4 

A7 

A6 

A5 


X = (A1 + 2* A2 + A3) - (A7 + 2* A6 + A5) 
Y = (A1 + 2* A8 + A7) - (A3 + 2* A4 + A5) 
G(X,Y) =\/ X 2 + Y 2 
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Based on tests with many nun-wave images, it was determined 
that the Sobel operator gave the best results for small objects. 
Examples using the edge detection algorithm are shown in Figs. 
3-2d through 3-5d , 3-8d, 3-9d and 3-10c. In the cloudy freeway 
scenes, Figs. 3 -2d and 3 -3d, this processing highlights the small 
vehicles and the edges of the roads. In Fig. 3-8d, the military 
vehicles stand out against the smooth background, making it easy 
to count them and to identify their formation. In the Long Beach 
Navy dockyard scene in Fig. 3-9, the image is very complicated, 
and the edge detection in this form does not aid in the analysis. 

One additional edge detection algorithm which is useful for 
distributed targets is the Marr-Hildreth technique (Ref. 10) 
which can be used to outline targets. This could be helpful with 
computers to automate the image analysis using object recognition 
algorithms. The first step in this technique is to smooth the 
image by convolving it with a Gaussian distribution: 

G(x, y) = exp(-(x 2 + y 2 )/2*Pi * a 2 ), (7) 

where a is the size of the smoothing area. The second step of 
this process is a convolution with the Laplacian of the two- 
dimensional Gaussian: 

2G(x, y) = -1/ (Pi*o 4 ) * (1 - (x 2 + y 2 ) *exp(- (x 2 + y 2 )/2cr 2 ) (8) 

The third step is to find the zero-crossings and draw a 
line. An example of this is shown with the ship image in Fig. 3- 
lOd. In this case the major structures of the ship are outlined, 
which could help in its identification. The limitation on this 
technique though, is that the targets must be large enough (~100 
pixels) to give recognizable shapes. A further disadvantage is 
that this operator is isotropic because of the symmetry of the 
Laplacian. Better results can be obtained with so-called 
anisotropic operators (for example see Ref. 11) . 
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Figure 3-9. Visual and rmi-wave images of the Long Beach na 
and the mm-wave data (b) cover an area of 0.5 
image is in (c) and the edge detected image in 
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Figure 3-10. MM-wave images of a large ship in the Long Beach, Ca. harbor just south of Terminal Island. 

The original data is in (a). The mm-wave data processed with the Laplace operator is in (b), 
the Sobel operator in (c) and the Marr-Hildreth algorithm in (d). 




CHAPTER IV. CONCLUSIONS 

In this study, techniques of color scale choice and 
adjustment on the brightness temperature scale, modification 
of the brightness temperature histogram, and edge detection 
techniques have been found to be very effective in the visual 
interpretation of the mm-wave images. Other techniques of 
smoothing and outlining were found not to aid in the analysis 
except in specialized cases. The display of the original mm- 
wave data along with the contrast enhanced and the edge 
detection images appears to be a good combination to assist 
an observer in the interpretation of this data. The observer 
can also adjust the color scale and gain factors using his 
intuition to bring out the details to understand the image. 
Having this information in real time is important in the 
usefulness of this data. Enhancing this data may lead to 
future steps to get computer recognition of the images. 
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APPENDIX 


Summary of the MM-Wave Image Processing Program Package 


To use this program package as summarized in the Table 
A-l, the following comments are necessary: 

- Files 

Each image consists of a header file - specified by xxxxxx 
INF - and a data file - specified by xxxxxx DAT. 

- Data File 

The image data are written with the JPL-subroutines DKREAD 
and DKWRIT (using corresponding VAX-VMS sub-routines) on and 
from disk. 

- Header File 

Each header consists of the following three numbers: number 

of image lines, number of samples per image line, and the 
specification of the data type (1* two byte integer data, 2* 
byte data) . The header data are written in the following 
FORTRAN - format: 316. 

- Greyman, Filter, Edge 

There are listed more functions in the menus than described 
in Chapter 3 . The mathematical background for these 
algorithms is described for Greyman in [4], for FILTER and 
EDGE in [ 3 ] . 
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Table A-l. Survey of Image Processing Programs 


FUncticn 

Display an Image 
loaded cn Disk 


load, Change and 
Manipulate Oder 
Table 

Program 





Name 

DISPIMAGE 


CSCDISP 


Man 






DISPLAY RAW IMAGES 

[1] 

EXIT 

[0] 


ZOCM OF IMAGES 

[2] 

ADJUST CENUR 



DISPLAY SUB-IMAGES 

[3] 

AND <3T CDICR 



DISPLAY A NEW IMA3I 

[4] 

M«P FRCM DISK 



SAVE 1034*1024 UPGES 

[5] 

EELE 

[2] 


DISPLAY 1024*1024 IMAGES 

[6] 

CKEAIE A NEW 



SAVE 1280*1024 DW2S 

[7] 

GDIGR MAP 

[3] 


DISPLAY 1280*1024 IMA2S 

[8] 

HIT CDICRMftP 



SAVE SUB-IMAO: FRCM 

[9] 

CN DISK 

[4] 


SCREEN CURSOR 

[10] 




exit 

[11] 



Rjncticn 






Apply Various Ocntrast 

Apply Various Filter Apply Various 


Enhanoanent tfethods to 

Methods to an Image 

Qiianoanerit Methods to 


an Image loaded cn Disk 

loaded cn Disk 

an Image loaded cn Disk 

Pruyidiu 





Name 



FEDUR 

edge; 


Maru 


***NCNUNE?R EUMJCEMENT*** 







SCEEL 

[1] 

HCBABILITY CISIRIHJriCN 

[1] 

3E0JHING 

[1] LAPLACE 

[2] 

EQUALIZATION 

[2] 

MEDIAN 

[2] RCBERIS 

[3] 

HYffRBCUZATICN 

[3] 

I0WEASS FI HIM* [3] PROCESS NWE DBGE [4] 

EXKNENIIAIIZATICN 

[4] 

HIGHEASS FEIIER[4] EXIT 

[5] 

RA YT FT CH-Fl NOTION 

[5] 

EANDfftSS FITIER[5] 


POIINZ-FLNCnCN 

[6] 

FRCCE3S NEW 





IMAGE 

[6] 




EXIT 

[7] 


***LENEAR ENHANCEMENT*** 





MIN-MAX HISICO^M SUEKH[7] 




5-85 ^HISTOGRAM SIFEICH 

[8] 




VARIABLE CUPPING 

[9] 




HCCESS NEW FILE 

[10] 




EXIT 

[11] 





